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ABSTRACT 

We present the 3D real space clustering power spectrum of a sample of ~ 600,000 
luminous red galaxies (LRGs) measured by the Sloan Digital Sky Survey (SDSS), using 
photometric redshifts. These galaxies are old, elliptical systems with strong 4000 A 
breaks, and have accurate photometric redshifts with an average error of Az — 0.03. 
This sample of galaxies ranges from redshift z = 0.2 to 0.6 over 3, 528 deg 2 of the sky, 
probing a volume of 1.5/i~ 3 Gpc 3 , making it the largest volume ever used for galaxy 
clustering measurements. We measure the angular clustering power spectrum in eight 
redshift slices and use well-calibrated redshift distributions to combine these into a high 
precision 3D real space power spectrum from k = 0.005/iMpc" 1 to k = 1/iMpc -1 . We 
detect power on gigaparsec scales, beyond the turnover in the matter power spectrum, 
at a ~ 2a significance for k < 0.01/iMpc -1 , increasing to 5.5a for k < 0.02/iMpc -1 . 
This detection of power is on scales significantly larger than those accessible to current 
spectroscopic redshift surveys. We also find evidence for baryonic oscillations, both in 
the power spectrum, as well as in fits to the baryon density, at a 2.5a confidence 
level. The large volume and resulting small statistical errors on the power spectrum 
allow us to constrain both the amplitude and scale dependence of the galaxy bias in 
cosmological fits. The statistical power of these data to constrain cosmology is ~ 1.7 
times better than previous clustering analyses. Varying the matter density and baryon 
fraction, we find Q M = 0.30 ± 0.03, and n b /n M = 0.18 ± 0.04, for a fixed Hubble 
constant of 70km/s/Mpc and a scale-invariant spectrum of initial perturbations. The 
detection of baryonic oscillations also allows us to measure the comoving distance 
to z = 0.5; we find a best fit distance of 1.73 ± 0.12Gpc, corresponding to a 6.5% 
error on the distance. These results demonstrate the ability to make precise clustering 
measurements with photometric surveys. 

1 INTRODUCTION 

The three dimensional distribution of galaxies has lone 
been recognized as a powerful cosmological probe jPeebled 
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density to have a simple relationship to the underlying mat- 
ter density; therefore, the clustering of galaxies is related 
to the clustering of the underlying matter. The two point 
correlation function of matter (or its Fourier transform, the 
power spectrum) is a sensitive probe of both the initial con- 
ditions of the Universe and its subsequent evolution. Indeed, 
if the matter density is well described by a Gaussian ran- 
dom field, then the power spectrum encodes all the informa- 
tion present in the field. It is therefore not surprising that 
a large fraction of the effort in observational cosmology has 
been devoted to measuring the spatial distribution of galax- 
ies, culminating in recent results fr om the Two-Deg ree Field 
Galaxy Redshift Survey (2dFGR S.ICole et all2005h an d the 
Sloan Digital Sky Survey (SDSS. lTegmark et alJl20oil . 

The spatial distribution of galaxies is also a standard 
ruler for cosmography. The expansion rate of the Universe 
as a function of redshift is a sensitive probe of its energy con- 
tent, and in particular, can be used to constrain the prop- 
erties of the "dark energy" resp onsible for the recent acce l- 
eration in the expansion (see eg. 151120051 : lEisensteinl feoOS ) . 
One approach to measure the expansion rate is to observe 
the apparent size of a standard ruler (and therefore, the an- 
gular diameter distance) at different redshifts to constrain 
the scale factor as a function of time. The power spectrum 
of the galaxy distribution has two features useful as stan- 
dard rulers. At k ~ 0.01/iMpc -1 , the power spectrum turns 
over from a k 1 slope (for a scale invariant spectrum of ini- 
tial fluctuations), to a k~ 3 spectrum, caused by modes that 
entered the horizon during radiation domination and were 
therefore suppressed. The precise position of this turnover 
is determined by the size of the horizon at matter-radiation 
equality, and corresponds to a physical scale determined 
by the matter (ftMh 2 ) and radiation densities (fijj/i 2 ). The 
other distinguishing feature is oscillations in the power spec- 
trum caused by acoustic waves in the baryon- photon plasma 
before hydrogen recombination at z ~ 1000 l|Peebles fc YvJ 

[ 197d: ISunvaev fc Zeldovichlll980t iBond fc Efstathioulll98l 
Holtzmanll989l : lEisenstein fc Hull998tlMeiksin et alll99sT . 
The physics of these oscillations are analogous to those of 
the cosmic microwave background, although their ampli- 
tude is suppressed because only ~ 1/6 of the matter in 
the Universe is composed of baryons. The scale of this fea- 
ture, again determined by the matter and radiation den- 
sities, is set by the sound horizon at hydrogen recombina- 
tion. This feature was first observed in ea rly 2005 both in 
the SPSS Lumino us Red Galaxy sample { Eisenstein e^jj 
l2005bl:lHutsl2006D and the 2dFGRS data ICole et all200Eft . 
Measuring the apparent size of both of these features at 
different redshifts opens up the possibility of directly mea- 
surin g the angular diameter distance as a function of red- 
shift (Eisenstein et al. 1998: ISeo fc Eisenstein| [2003: Linder 
20031: llVIatsubara fc Szalavll2003l iBlake fc Glazebrookll2003t 
Hu fc Haimanl 120031: iMatsubaral [20041: ISeo fc Eisenstein! 
20051: IWhitd 120051: IBlake fc Bridldl^Ool IBlake et alJl2006t 
Dolnev et alJl2006l) . 

Traditionally, measurements of galaxy clustering rely 
on spectroscopic redshifts to estimate distances to galax- 
ies. Even with modern CCDs and high throughput multi- 



fiber spectrographs, acquiring them is an expensive, time- 
consuming process compared with just imaging the sky. For 
instance, the SDSS spends about one-fifth of the time imag- 
ing the sky, and the rest on spectroscopy. Furthermore, the 
ultimate accuracy of distance estimates from spectroscopy 
is limited by peculiar velocities of ~ lOOOkm/s, a signifi- 
cant mismatch with the intrinsic spectroscopic accuracy of 
~ lOkm/s. 

Large multi-band imaging surveys allow for the possi- 
bility of replacing spectroscopic with photometric redshifts. 
The advantage is relative efficiency of imaging over spec- 
troscopy. Given a constant amount of telescope time, one 
can image both wider areas and deeper volumes than would 
be possible with spectroscopy, allowing one to probe both 
larger scales and larger volumes. Fu rthermore, the accu- 
racy o f photometric distance estimates llPadmanabhan et alJ 
l2005al) . cAz ~ 10, OOOkm/s is more closely matched (al- 
though still not optimal) to the intrinisic uncertainties in 
the distance-redshift relations. 

One aim of this paper is to demonstrate the practicality 
of such an approach by applying it to real data. We start 
with the five band imaging of the SDSS, and photometri- 
cally select a sample of luminous red galaxies; these galaxies 
have a strong 4000 A break in their spectral energy distri- 
butions, making uniform selection and accurate photometric 
redshifts possible. We then measure the angular clustering 
power spectrum as a function of redshift, and "stack" these 
individual 2D power spectra to obtain an estimate of the 3D 
clustering power spectrum. Using the photometric survey al- 
lows us to probe both larger scales and higher redshifts than 
is possible with the SDSS spectroscopic samples. 

We pay special attention to the systematics unique 
to photometric surveys, and develop techniques to test for 
these. Stellar contamination, variations in star-galaxy sepa- 
ration with seeing, uncertainties in Galactic extinction, and 
variations in the photometric calibration all can masquerade 
as large scale structure, making it essential to understand 
the extent of their contamination. Furthermore, stacking the 
angular power spectra to measure the 3D clustering of galax- 
ies requires testing our understanding of the photometric 
redshifts and their errors. 

The paper is organized as follows : Sec. ^describes the 
construction of the sample; Sec. [3] then discusses the mea- 
surement of the angular power spectrum and the associ- 
ated checks for systematics. These angular power spectra are 
then stacked to estimate the 3D power spectrum (Sec. [IJ, 
and preliminary cosmological parameters are estimated in 
Sec.|S] We conclude in Sec.|H] Wherever not explicitly men- 
tioned, we assume a flat ACDM cosmology with Qm = 0.3, 
ilb = 0.05, h — 0.7, a scale invariant primordial power spec- 
trum, and as = 0.9. 

2 THE SAMPLE 
2.1 The Data 

The Sloan Digital Sky Survey J York et alJl2000h is an on- 
going effort to image approximately n steradians of the 
sky, and obtain spec tra of approximately one million of 
the d etected objects iStrauss et al.l 120021: | Eisen stein et alJ 
1200 ll) . The imaging is carried out by drift-scanning the 
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sky in photomet ric condition s llHogg et all 1200 lli . using a 
2.5m telescope jG^mi^^ilJj2006 ) in fi ve bands (ugriz) 
jFukueita et al]ll996l : ISmith et alll2002Tl usin g a specially 
designed wide-field came7a ^l&m^ie^nTll99'^. Using these 
data, objects are targeted for spectroscopy ^Richards et alJ 
120021: iBlanton et"aD 120031) and are observed with a 640- 
fiber spectrograph on the same telescope. All of these 
data are processed by completely automated pipelines that 
detect and measure photometric pr operties of objects, 
and astrometrically calibrate the d ata llLupton et aljl200lt 
iPier et aljEocB llvezic et al.ll2004l) . The first phase of the 
SDSS is com plete and has produced five major data re- 
leases dStouehton et alJ 120021: lAbazaiian et alJl20o3. 12004L 
lAdel 



120051: lAdelman- McCarthy et al.ll2006ft 1 . This paper uses all 
data observed through Fall 2003 (corresponding approxi- 
mately to SDSS Data Release 3), reduced as described by 
iFinkbeiner et alJ 12004^1 . 



2.2 Photometric Calibration 

Measurements of large scale structure with a photometric 
survey require uniform photometric calibrations over the en- 
tire survey region. Traditional methods of calibrating imag- 
ing data involve comparisons with secondary "standard" 
stars. The precision of such methods is limited by transfor- 
mations between different photometric systems, and there 
is no control over the relative photometry over the entire 
survey region. The approach we adopt with these data is 
to use repeat observations of stars to constrain the pho- 
tometric calibration of SDSS " runs" , analogou s to CMB 
map-making techniques (see eg. lTegmarkl ll997aV Since all 
observations are made with the same telescope, there are 
none of the uncertainties associated with using auxiliary 
data. Also, using overlaps allows one to control the rela- 
tive calibration over connected regions of survey. The only 
uncertainty is the overall zeropoint of the survey, which we 
match to published SDSS cal ibrations. The ab o ve me thod 
has been briefly des cribed by IFinkbeiner et all <l2004h and 
IBlanton et alJ J2005T) . and will be explained in detail in a 
future publication. 



2.3 Denning Luminous Red Galaxies 

Tracers of the large scale structure of the Universe must sat- 
isfy a number of criteria. They must probe a large cosmolog- 
ical volume to overcome sample variance, and have a high 
number density so shot noise is sub-dominant on the scales 
of interest. Furthermore, it must be possible to uniformly 
select these galaxies over the entire volume of interest. Fi- 
nally, if spectroscopic redshifts are unavailable, they should 
have well characterized photometric redshifts (and errors), 
and redshift distributions. 

The usefulness of LR Gs as a cosmologi c al probe has long 
been appreciated (jGladders fc Yed l2000i lEiscnstci n et alJ 
1200 J) . These are typically the most luminous galaxies in 
the Universe, and therefore probe cosmologically interest- 
ing volumes. In addition, these galaxies are generically 
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Figure 1. A model spect rum of an elliptical galaxy, taken from 
Bruzual & Chariot (2003), shown at three redshifts. The model 
assumes a single burst of star formation 11 Gyr ago and solar 
metallicity; the effect of evolution is not shown for simplicity. Also 
overplotted are the response functions (including atmospheric ab- 
sorption) for the five SDSS filters. The break in the spectrum at 
4000 A, and its migration through the SDSS filters is clearly seen. 



old stellar systems with uniform spectral energy distri- 
butions (SEDs) characterized principally by a strong dis- 
continuity at 4000 A (Fig. 0. This combination of a uni- 
form SED and a strong 4000 A break make LRGs an 
ideal candidate for photometric reds hift algorithms, with 
redshi ft accuracies of a z ~ 0.03 ijPadmanabhan et alJ 
2005a ). LRGs have been used for a number of studies 
torata et alJ 120041: lEisenstein et alJ l2005al IZehavi et alJ 
2005: IPadmanabhan et al.ll2005blh ' most notably for the de- 
tection of the bar yonic acoustic peak in th e galaxy autocor- 
relation function feisenstein et al1l2005tJ) 

The photometri c selection criteria w e adopt were dis- 
cussed in detail by IPadmanabhan et, alJ (|2005ar ) and are 
summarized below. We start with a model spectrum of an 
early type gala xy from the stellar population synthesis mod- 
els of Ifjruzual fc CharlotJ feOOSt (Fig. Q. This particular 
spectrum is derived from a single burst of star formation 
11 Gyr ago (implying a redshift of formation, Zf orm ~ 2.6), 
evolved to the present, and is typical of LRG spectra. In 
particular, the 4000 A break is very prominent. To motivate 
our selection criteria, we passively evolve this spectrum in 
redshift (taking the evolution of the strength of the 4000 
A break into account), and project it through the SDSS 
filters; the resulting colour track in g — r — i space as a func- 
tion of redshift is shown in Fig. |5] The bend in the track 
around z ~ 0.4, as the 4000 A break redshifts from the g 
to r band, naturally suggests two selection criteria - a low 
redshift sample (Cut I), nominally from z ~ 0.2 — 0.4, and 
a high redshift sample (Cut II), from z ~ .4 — 0.6. We 
define the two colours feisenstein et all 1200 ll and private 
commun.) 
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Figure 2. The top panel shows simulated g — r and r — i colours 
of an early-type galaxy as a function of redshift. The spectrum 
used to generate the track is the same as in Fig. but evolved 
in redshift. Also shown are the colour cuts for Cut I (dashed, 
black) and Cut II galaxies (solid, blue). The lower panel shows 
the colours en (diamonds, black) and (triangles, red), as a 
function of redshift. Also shown are fiducial redshift boundaries 
for Cut I (0.2 - 0.4) and Cut II (0.4 - 0.6). Note that the range 
in g — r is identical to the range in 1 + z. 



d x = {r-i)-{g-r)/%^r-i . (2) 

We now make the following colour selections, 

Cut I : | c± |< 0.2 ; (3) 

Cut II : d x > 0.55 , (4) 

g - r > 1.4 , (5) 



as shown in Fig. |21 The final cut, g — r > 1.4, isolates our 
sample from the stellar locus. In addition to these selection 
criteria, we eliminate all galaxies with g — r > 3 and r — 
i > 1.5; these constraints eliminate no real galaxies, but are 

effective at removing stars with unusual co lours. 

Unfortunately, as emphasized in lEisenstein et alJ 

i200ll) . these simple colour cuts are not sufficient to select 
LRGs due to an accidental degeneracy in the SDSS filters 
that causes all galaxies, irrespective of type, to lie very close 
to the low redshift early type locus. We therefore follow the 
discussion there and impose a cut in absolute magnitude. 
We implement this by defining a colour as a proxy for red- 
shift and then translating the absolute magnitude cut into a 
colour-apparent magnitude cut. We see from Fig. |2]that d± 
correlates strongly with redshift and is appropriate to use 
for Cut II. For Cut I, we define, 

c N = 0.7(5 - r) + 1.2(r - i - 0.18) , (6) 

which is approximately parallel to the low redshift locus. 
Given these, we further impose 

Cut I : rp e tro < 13.6 +C| |/0.3 , 

rpetro < 19.7 ; (7) 




Figure 3. The angular selection function of the LRGs with 
the "Northern Celestial Cap" (black) and the "Equatorial Cap" 
(blue) shown. The lightly shaded (green) region of the Equato- 
rial cap (b < 45° , shown as a dashed line) is excluded because of 
possible stellar contamination. The gaps in the selection function 
are due to missing data and exclusion around bright stars. Also 
shown is the Galactic equator (solid line). 

Cut II : i < 18.3 + 2d ± , 

i < 20 . (8) 

Note we use the r band Petrosian magnitude (rp et ro) for 
consistency with the SDSS LRG target selection. We note 
that Cut I is identical (except for the magnitude cuts in Eqs. 
to the SDSS LRG Cut I, while Cut II was chosen to yield a 
population consistent with Cut I (see below). This was inten- 
tionally done to maximize the overlap between any sample 
selected using these cuts, and the SDSS LRG spectroscopic 
sample. The switch to the i band for Cut II also requires 
explanation. As is clear from FigQ the 4000 A break is red- 
shifting through the r band throughout the fiducial redshift 
range of Cut II. This implies that the K-corrections to the r 
band are very sensitive to redshift; the i band K-corrections 
are much less sensitive to redshift allowing for a more robust 
selection. 

Finally, we augment the star-galaxy separation from 
SDSS with the following cuts designed to minimize stellar 
contamination, 

Cut I : r PSF ~r> 0.3, 
Cut II: i P SF - i > 0.2(21 - i) , 

r de v > 0.2 , (9) 

where (r,i)psF are the SDSS PSF magnitudes, while rdev 
is the deVaucouleurs radius of the galaxy in arcseconds. 

2.4 Angular and Redshift Distributions 

Applying the above selection criteria to the ~ 5500 degrees 
of photometric SDSS imaging considered in this paper yields 
a catalog of approximately 900,000 galaxies. We pixelize 
these galaxies as a num ber overdensity, 5 a = Sn/n, onto 
a HEALPix pixelization dGorski et alj ll999) of the sphere, 
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Figure 4. (Top) The photometric redshift distribution of the 
LRG sample. (Bottom) The deconvolved selection functions for 
photometric redshift slices with Az = 0.05 from z = 0.2 to z = 
0.6. The dotted lines are the mean redshifts of each of the slices. 

with 3,145,728 pixels (HEALPix resolution 9). We exclude 
regio ns where the extinction in the r-band iSchleeel et alJ 
1998) exceeds 0.2 magnitudes, as well as maski ng regions 
aroun d stars in the Tycho astrometric catalog ijHog et alJ 
2000). We also exclude data from the three southern SDSS 
stripes due to difficulties in photometrically calibrating them 
relative to the rest of the data, due to the lack of any overlap. 
The resulting angular selection function is shown in Fig. |3] 
The angular coverage naturally divides into two regions, 
which we refer to as the "Northern Celestial Cap" (NCC) 
and the "Equatorial Cap" (EC), based on their positions 
on the celestial sphere. As discussed below, we additionally 
excise regions in the EC with b < 45° due to possible stel- 
lar contamination. The final angular selection function cov- 
ers a solid angle of 2,384 square degrees (181,766 resolution 
9 HEALPix pixels) in the NCC, and 1,144 square degrees 
(87,263 resolution 9 HEALPix pixels) in the EC. 

The calibration and accuracy of photometric redshift 
algorithm s for this sample h ave been discussed in detail by 
iPadmanabhan et alJ i2005al) . We compute photometric red- 
shifts for all the galaxies in the sample using the simple 
template fitting algorithm described there; these redshifts 
have calibrated errors of a z ~ 0.025 at z ~ 0.2 that in- 
crease to g z ~ 0.05 at z ~ 0.6. The resulting photometric 
redshift distribution is in Fig. 2] The sample is divided into 
8 photometric redshift slices of thickness Az = 0.05 (zOO 
through 207), and the underlying redshift distributions for 
each slice are estimated using the deconvolution algorithm 
presented in the above reference. These redshift distribu- 
tions are plotted in Fig. 2] while properties of the different 
slices are summarized in Table Q 

2.5 Sample Systematics 

There are a number of systematic effects in photometric 
samples that contaminate clustering - stellar contamination, 
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0.275 


0.276 


20377 
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1.52± 0.06 


z02 


0.325 


0.326 


21759 


10768 


1.67± 0.07 


z03 


0.375 


0.376 


28345 


12706 


1.94± 0.06 


z04 


0.425 


0.445 


41527 


18767 


1.75± 0.06 


z05 


0.475 


0.506 


71131 


33000 


1.73± 0.04 


z06 


0.525 


0.552 


65324 


30281 


1.80± 0.04 


z()7 


0.575 


0.602 


46185 


20504 


1.85± 0.05 



Table 1. Descriptions of the 8 Az = 0.05 redshift slices; z mic i 
is the midpoint of the redshift interval, while z m ean is the mean 
redshift of the slice. Also listed are the number of galaxies (N ga i) 
for the "Northern Celestial Cap" (NCC), and the "Equatorial 
Cap" (EC), and the linear bias of each redshift slice, b g . 
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Figure 5. The average number of LRGs per resolution 6 
HEALPix pixel (approximately 1 deg 2 in area) as a function of 
Galactic latitude, for the two disjoint caps. The contours are 16%, 
50% and 84%. There is some evidence for stellar contamination 
(sec text for more details) at low Galactic latitudes in the Equa- 
torial Cap; excising the region b < 45° removes the problematic 
regions. 



angular and radial modulation of the selection due to see- 
ing variations, extinction, and errors in our modelling of the 
galaxy population. Fig. plots the areal LRG density as a 
function of Galactic latitude; one would expect any leakage 
in the star-galaxy separation to increase at lower latitudes 
where the stellar density is higher. We see no increase for 
the NCC, but observe an increase for b < 45° for the EC. 
This is further borne out by Fig. |S| where we plot the LRG 
density versus the density of stars with SDSS PSF magni- 
tudes 18.0 < rpsF < 19.5, where the magnitude limits were 
chosen so that the SDSS star-galaxy separation is essentially 
perfect. Although the precise signature of such contamina- 
tion on the clustering signal is unclear, we choose to be con- 
servative and exclude regions below b — 45° (Figs^and 
this reduces the area of the EC by 25%. 
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Figure 6. The same as Fig.[S]except now as a function of stellar 
density. The stellar density is estimated analogous to the galaxy 
density, selecting stars with PSF fluxes rpgp between 18.0 and 
19.5. The vertical line shows the position of the b < 45° cut in 
the Equatorial Cap. 

To understand the nature of this contamination, we 
consider the subset of galaxies for which SDSS has mea- 
sured spectra. We find that 118,053 (13.1%) galaxies in 
the photometric sample have measured spectra. Of these, 
662 (0.56%) are unambiguously classified as stars (475 ob- 
jects) or quasars (187 objects). The quasars are at low 
(0.1 < z < 0.25) redshifts, while the stars are almost en- 
tirely K and M stars, and are preferentially at lower Galactic 
latitudes, consistent with the above. Inspecting the imaging 
data shows that these are either late-type stars blended with 
other stars (approximately 2/3), late-type stars blended 
with background galaxies (approximately 1/3), and a smat- 
tering of star-artefact blends. Note that this explains the 
dependence with Galactic latitude and stellar density; one 
would naively expect the number of star-star blends to scale 
as the square of the stellar density, while the star-galaxy 
and star-artefact blends should roughly scale as the stel- 
lar density. We emphasize that the levels of contamination 
obtained this way are approximate, since the spectroscopic 
survey has a brighter apparent luminosity limit than our 
photometric catalog, and the contamination could increase 
with decreasing luminosity. 

To test for the possible modulation of the LRG selec- 
tion due to angular variations in seeing and extinction, we 
consider the areal density of LRGs observed as a function 
of seeing (as measured by the FW HM of the r-band PSF) 
and extinction (Schlcgcl et al. 1998). These distributions are 
plotted in Figs.[7|and |H| We find that the density is constant 
to ~ 2% over most of the range of seeing and extinction in 
the survey. We do observe deviations at the very edges of 
the distributions, but the total area with these extremes in 
seeing and extinction is negligible (as seen in the top panels 
of the figures), and therefore, do not affect clustering mea- 
surements. 

Finally, to test sample uniformity as a function of red- 



Figure 7. (top) The histogram is the (normalized) distribution 
of galaxies as a function of the PSF FWHM (measured in arcsec- 
onds) in the r band. The (red) curve is the fraction of the total 
survey area with the same PSF FWHM. The agreement between 
them suggests that the galaxy selection algorithm is unaffected 
by seeing, (bottom) The galaxy surface density as a function of 
seeing. The two distributions are identical at the 2% level except 
at the edges where the relevant survey area is negligible. 
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Figure 8. Analogous to Fig. 13 except with extinction in the r- 
band from Schlegel et al. 119981) . We truncate at A r = 0.2 corre- 
sponding to the cut in the angular selection function. 

shift, we consider the luminosity distribution as a function of 
redshift. A constant luminosity distribution over the redshift 
range would suggest that we were selecting comparable pop- 
ulations of galaxies. A complication is that we must use pho- 
tometric redshifts to compute absolute magnitudes; biases 
in the photometric redshifts could alter the inferred luminos- 
ity distributions. We estimate the magnitude of such biases 
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-23.0 




Figure 9. The conditional i— band magnitude distribution as a 
function of redshift. The absolute magnitude is computed assum- 
ing the photometric redshift. The contours show the 16%, 50%, 
and 84% levels, while the horizontal (red) line is the median mag- 
nitude for the entire sample. The glitch at z ~ 0.4 corresponds to 
the transition between Cut I and Cut II LRGs, while the increase 
at z ~ 0.55 is due to the magnitude limit in Cut II. 



from Table At low redshifts, the photometric redshifts 
are essentially unbiased, whereas at high redshifts, the pho- 
tometric redshifts underestimate the true redshift by about 
Az = 0.025, which translates into an overestimation of the 
magnitude by about AM = 0.1 — 0.15 magnitudes. 

The observed conditional luminosity distribution as a 
function of redshift is in Fig. |2l The median luminosity is 
constant to approximately AM = 0.1 over the redshift range 
of interest, with a width of ~ 0.7 magnitudes (compared 
with a potential bias of 0.15 magnitudes above). The distri- 
bution has two distinguishing features, a glitch at z ~ 0.4 
and increasing luminosities at higher redshifts. The glitch 
at z ~ 0.4 corresponds to the transition between Cut I 
and Cut II at the point where colour tracks bend sharply 
in Fig. |21 and highlights a difficulty in uniformly select- 
ing galaxies in that region. The increase in luminosities at 
~ 2 = 0.55 is due to the magnitude limits imposed in Cut II. 
Except for these two features, we conclude that our selection 
criteria yield an approximately uniform galaxy population 
from z = 0.2 to z — 0.55. 



3 THE ANGULAR POWER SPECTRUM 
3.1 Projections on the sky 

We relate the projected angular power spectrum to the 
underlying three dimensiona l power spectrum; o ur deriva- 
tion follows the discu ssion in lHuterer et alJ J200ll) (see also 
iTeemark et al.l I2002L and references therein). We describe 
the galaxy distribution by an isotropic 3D density field, 
S g ,3D, and its power spectrum P(k) denned by, 



(10) 



Projecting this density field on the sky along n, we obtain 
1 



J dy <j>{y) 



dy<t>(y)S g ,3D(y,yn) 



(11) 



where y is the comoving distance, and 4>{y) is the radial se- 
lection function. For now, we ignore the effect of peculiar 
velocities, and therefore do not distinguish between real and 
redshift space quantities. Fourier transforming the 3D den- 
sity field and making use of the identity, 



(12) 



we obtain, 



S g (h) = / dyf(y) 



d*k 



(2tt) 3 

OO 

x^i i (2/ + l)j i (fcy)P ; (n-k) 



<5 9 ,3D(y,k) 



(13) 



where ji (x) and Pi (x) are the i* order spherical Bessel func- 
tions and Legendre polynomials respectively. We define the 
weighting function, f(y) by 

4>(y) 



J dycj>{y) 



(14) 



Since the density field is isotropic, we expand it in Legendre 
polynomials to obtain, 



5 g ,i =i l dyf(y) 



d s k 
(2tt) 3 



8 a ^ D {y,)^)ji{ky) 



(15) 



In order to proceed, we assume that the selection func- 
tion is narrow in redshift, allowing us to ignore the evolution 
of the density field. The above equation can then be written 

as 



8g,i = i l 



d A k 
(2tt) 3 



* B ,sx>(k)Wi(A) , 



(16) 



where we implicitly assume that the density field is at the 
median redshift of the selection function. The window func- 
tion, Wi(k), describes the mapping of k to / and is given 
by, 



Wi(k)= / dyf(y)ji(ky) 



(17) 



It is now straightforward to compute the angular power spec- 
trum, 



( : = 'A,,<\, ; / ^-A 2 (k)W?(k) , 



(18) 



where A 2 (fe) is the variance per logarithmic wavenumber, 



A\k) 



-Ank 3 P(k) 



(19) 



(2tt) 3 

Similarly, the cross correlation between two selection func- 
tions, cpi and 02, is give by 



(20) 



We have not distinguished between the galaxy and mat- 
ter power spectrum above. On large scales, we simply assume 

Pg(k) = b 2 g P(k) , (21) 
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Figure 10. The theoretical angular power spectra for each of the 
redshift slices in Fig.|l] The heavy and light solid lines show non- 
linear and linear auto power spectra, while the dashed lines show 
the cross power spectra with the adjacent slice at higher redshift. 
The dotted lines show the effect of redshift space distortions on 
both the auto and cross power spectra, assuming /3 = 0.3. The 
panel z04 also shows the cross correlation power with z06 and 
z07. 



where P g {k) and P(k) are the galaxy and matter power spec- 
tra respectively, and b g is the linear galaxy bias. This is a 
good approximation on large scales ijScherrer fc Weinberg! 
1998), but breaks down on smaller scales; we defer a dis- 
cussion of its regime of validity, as well as the nonlinear 
evolution of the power spectrum to a later section. 

Fig. 1101 shows the predicted angular power spectra for 
the eight redshift distributions in Fig. 2] assuming our fidu- 
cial cosmology; also shown are the cross-correlation power 
spectra for adjacent s lices. We assume b a = 1, and use the 
halof it prescription jSmith et alJl2003F) to evolve the mat- 
ter power spectrum into the nonlinear regime. The increase 
in the amplitude of the power spectrum on large scales (low 
L) with decreasing redshift is due to the linear growth of 
structure, while the increase in power on small scales (large 
L) is from the nonlinear collapse of structures. The "knee" in 
the power spectrum between L ~ 10 — 30 corresponds to the 
turnover in the 3D power spectrum P(k), where the shape 
changes from P(k) ~ k to P(k) ~ k~ 3 (in the linear regime). 
This scale corresponds to the horizon at matter-radiation 
equality and is constant in comoving coordinates. However, 
with increasing radial distances to the redshift slices, the 
apparent angular size decreases with redshift, and we see 
the "knee" shift from low L (large angular scales) at low 
redshifts to high L (small angular scales) at high redshifts. 
This illustrates the potential use of the power spectrum as a 
standard ruler for cosmography; given the size of the horizon 



at matter-radiation equality (independent of dark energy), 
one can probe the evolution of the universe during the dark 
energy dominated phase. A second such standard ruler is 
the baryonic oscillations in the matter power spectrum visi- 
ble at L ~ 100. However, its amplitude is suppressed in the 
individual angular power spectra by the smoothing due to 
the thickness of the redshift slices. 

Finally, we note that the cross-correlation between ad- 
jacent slices is non-negligible. This is easily understood by 
considering Fig. 01 where we note that there is considerable 
overlap between adjacent slices. Furthermore, this overlap 
increases with increasing redshift due to larger photometric 
redshift errors; this too is reflected in the cross-correlations. 
Going to more widely separated slices reduces the cross- 
correlation due to smaller overlaps. Note that the level of 
correlation seen in Fig. 1101 is only true on large scales; on 
smaller scales, uncorrelated Poisson noise (since the galaxy 
samples are disjoint) erases these correlations. 



3.1.1 Redshift Space Distortions 

The above discussion ignored the effect of peculiar veloci- 
ties on the observed clustering power spectrum. For broad 
redshift selection functions, the projection on to the sphere 
erases redshift space distortions; however, as the selection 
function becomes narrow, they become more important. 
We calculate their effect below, following the formalism of 
iFisher et aD Jl994) . 

We start with Eq. [TTJ 



1 + S g (h)= / dyf(s)[l + 6 g , 3D (y,yn)] 



(22) 



where we have now written the weighting function as a func- 
tion of redshift distance, s = y + v • n, and we have left the 
monopole contribution to the projected galaxy density ex- 
plicit. Assuming the peculiar velocities are small compared 
with the thickness of the redshift slice, we Taylor expand 
the weight function to linear order, 



/00 «/(!/) + f (v(yn)-n) 
dy 



(23) 



Substituting this expression into Eq. 1221 we note that at 
linear order, redshift space distortions only imprint fluctua- 
tions on the monopole component of the galaxy density. This 
allows us to separate the 2D galaxy density into two terms, 
S g = 8g + Sg, where 5° is the term discussed above, while 8 r g 
are the redshift space distortions. Fourier transforming the 
velocity field, we find that, 



dy 



d J k 



dy J (2ttY 



v(k). 



-zk-ny 



(24) 



The linearized continuity equation allows us to relate the 
velocity and density perturbations, 



v(k) = 



(25) 



where j3 is the redshift distortion parameter given approxi- 
mately by f3 « n„ 6 /6 9 . Substituting Eq. 03 into Eq.l2lland 
taking the Legendre transform, we can rewrite this equation 
in the form of Eq. 1161 where the window function now has 
an additional component given by, 
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Wf(k) 



dy -^fi( k y) 



(26) 



where j[ is the derivative of the spherical Bessel function 
with respect to its argument. By a repeated application of 
the recurrence relation lji-i — (I + l)ji+i — (2' + l)j'i, and 
integrating by parts, 



W[(k)=(3 / dyf(y) 



(2l 2 + 21-1) 



1(1-1) 



(21 - l)(2l + 1) 



ji-2(ky) 



(2l + 3)(2l 

(l + l)(l + 2) 
(2l + l)(2l + 3) 



ji(ky) 



-1) 
ji+2(ky) 



(27) 



It is interesting to note that this result could have been 
equivalently derived by starting from the iKaiserl il987t) en- 
hancement of the 3D power spectrum due to redshift space 
distortions, P g (k) — > P s (k)[l + /3(k • n)], and integrating 
along the line of sight as in Sec. 13.11 the I ± 2 spherical 
Bessel functions result from the coupling of the k ■ n an- 
gular dependence to the Legendre polynomials. Also inter- 
esting is the I ^> limit of the above equation; for I suffi- 
ciently large, J dyf(y)ji(y) « J dy f(y)ji± 2 (y), and W r (k) 
vanishes. Physically, this corresponds to the radial velocity 
perturbations being erased by the projection on to the sky. 

Fig. 1101 shows the effects of redshift space distortions 
on the angular power spectra for the eight redshift slices we 
are considering. Note that they contribute significantly only 
on the largest scales (I <~ 30), justifying our use of linear 
theory. 



3.2 Power Spectrum Estimation 

The theory behind optimal power spectrum estimation is 
now well established, and so we limit ourselves to details spe- 
cific to this discussion, an d refer the reader to the numerous 
references on the subject jTeemark et al 1 ll99d : ISeliak1ll998t 
iPadmanabhan et al.ll20oH and references therein). 

We start by parametrizing the power spectrum with 
twenty step functions in I, CI, 



ViCl 



(28) 



where the p; are the parameters that determine the power 
spectrum. We form quadratic combinations of the data, 



C[i — —X CiC C^X 



(29) 



where x is a vector of pixelized galaxy overdensities, C is 
the covariance matrix of the data, and d is the derivative 
of the covariance matrix with respect to pi. The covariance 
matrix requires a prior power spectrum to account for cosmic 
variance; we estimate the prior by computing an estimate of 
the power spectrum with a flat prior and then iterating once. 
We also construct the Fisher matrix, 



F y = itr [C.C^C.C- 1 ] 



(30) 



covariance matrix F . 

A final note on implementation - the dimension of the 
data covariance matrix is given by the number of pixels 
in the data. This quickly makes any direct implementation 
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Figure 11. (top) The average recovered power spectrum from 
100 simulated realizations. The dashed (red) line is the input 
power spectrum. The error bars are the errors per realization, 
and are not the error on the mean. Note that we have suppressed 
the input power spectrum by a constant factor relative to the 
expected power to avoid getting 5 > 1; the galaxy number density 
was boosted by the same factor to reduce the shot noise, (bottom) 
X comparing the errors on the power spectrum derived from the 
Fisher matrix versus those obtained from the run to run variance 
of the simulations for each of the eight redshift slices and the two 
angular caps. Assuming the density field is Gaussian, the error on 
the power spectrum errors (a) is <r/\ / 2N, where N is the number 
of the simulations. Also shown are the ltr and 3cr lines. 



of this algorithm impractical. We the r efore use the algo- 
rithm outlined bv |Padmanabliaii_et_alJ j2003j) . modified for 
a spherical geometry as in iHirata et ail (120041) . 



3.3 Simulations 

Before applying the above algorithm to the LRG catalog, 
we apply it to simulated data. In addition to testing the ac- 
curacy of our power spectrum code, we would also like to 
understand the correlations between the NCC and the EC, 
allowing us to combine separate power spectrum measure- 
ments. 

In order to do so, we use the prior power spectra for 
each redshift slice to simulate a Gaussian random field over 
the entire sphere. We then Poisson distribute galaxies with 
probability (1 + S)/2 over the survey region, trimmed with 
the angular selection funct ion. One technical complication 
iPadmanabhan et al.ll200ll) is that the measured amplitude 
of the power spectrum results in a number of points with 
\S\ > 1, making simple Poisson sampling impossible. To 
avoid this, we suppress the power spectrum by a constant 
factor, and boost the number density of galaxies by the same 
factor to ensure that the shot noise is similarly suppressed. 
We generate 100 such simulations for the eight redshift slices 
and both angular caps separately, matching the observed 
numbers of galaxies in each case; although the different red- 
shift bins are uncorrelated, the angular caps are based on 



© 0000 RAS, MNRAS 000, 000-000 



10 Padmanabhan et al 



correlated density fields. This allows us to estimate the co- 
variance between power spectrum measurements made for 
the different caps. Our goal here is not to realistically simu- 
late galaxy formation, but to test our pipelines, and the re- 
sulting measurements and errors; Gaussian simulations are 
sufficient for this purpose. 

The results from one set of fOO simulations are shown 
in the top panel of Fig. If f I the recovered power spectrum 
agrees well with the input power spectrum. The bottom 
panel of the same plot compares the errors as measured by 
the inverse of the Fisher matrix with those obtained from 
the run to run variance of the simulations. Assuming Gaus- 
sianity, these errors should themselves have a relative error 
given by Aa/a — l/V2iV where N = fOO is the number of 
the simulations. As evident from the figure, the run to run 
variance agrees (within the expected errors) with the errors 
from the Fisher matrix. 



3.4 The Angular Power Spectrum 

Fig. ll2l shows the measured angular power spectrum for the 
eight redshift slices, with the two angular caps being mea- 
sured separately. The difficulty with processing the two an- 
gular caps simultaneously is that errors in photometric cali- 
bration masquerade as large scale power. While it is possible 
to control these systematics in regions with overlaps in the 
data, the two angular caps are disconnected; therefore, any 
relative calibration between the two caps must be indirect 
(eg. considering data taken on the same night, and assuming 
that the calibration is constant through the night). Unfor- 
tunately, the expected power on these scales is also small 
(A 2 ~ 10~ 3 ), and so we choose to be conservative and mea- 
sure the angular power spectrum for the caps separately. We 
combine these using the simulations of the previous section 
to correctly take the covariances between the two caps into 
account, fn order to avoid mixing power between different 
angular scales, we simply use constant weights proportional 
to the area (0.67 and 0.33 for the NCC and EC respec- 
tively); these are approximately the same weights that one 
would have obtained by inverse variance weighting. The final 
results are in Fig. If 21 

3.5 Bias and (3 

An immediate question is whether the power spectra in 
Fig. If 21 are consistent with being derived from a single 3D 
power spectrum, appropriately normalized to account for 
bias and the evolution of structure. We start with the linear 
3D power spectrum for our fiducial cosmology, and project 
it to a 2D power spectrum C^ gg , using the formalism of 
Sec. 13. II We also compute the effect of redshift space distor- 
tions, whose normalization we parametrize by f3, along the 
lines of Sec. 13.1.11 giving us two more power spectra, C^ gv 
and Ci yVV . The total power spectrum is, 

Ci = b 2 g (Ci, gg + 2pCi tBV + f3 2 Ci, vv ) , (31) 

where b g is the linear bias of the LRGs. The three power 
spectra represent correlations of the galaxy density with it- 
self (gg), the velocity perturbations (the source of linear red- 
shift distortions) with itself (vv), and their cross correlation 
(gv) . We also note (as emphasized in Sec. 13.1.11 that the 
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Figure 13. Bias as function of redshift, as estimated for the 
eight redshift slices, marginalizing over redshift space distortions. 
Note that the fourth slice, with its anomalous bias, corresponds 
in redshift to the glitch seen in Fig. |2] The dashed line shows the 
best linear fit to the all eight bias values, while the dotted line 
excludes the fourth data point. Note that we have ignored the 
correlations between the different redshift slices for the fit. 

redshift distortions only affect the largest scales; therefore, 
the linear assumption is justified. We can now explore the 
X 2 likelihood surface as a function of b and (3 for each of the 
redshift slices. In practice, j3 is not strongly constrained by 
these data, and so we marginalize over it when estimating 
the bias. 

The best fit models are compared with the data in 
Fig. 1121 while the bias for the eight redshift slices is in 
Fig. 1131 We do not fit to the entire power spectrum, but limit 
ourselves to scales larger than the nominal nonlinear cutoff 
at k = 0.1/i.Mpc , the angular scales corresponding to this 
restriction are marked in Fig. 1121 Our starting hypothesis 
- that the angular power spectra are derived from a single 
3D power spectrum - appears to be well motivated. Inter- 
estingly, the halofit nonlinear prescription for the matter 
power spectrum fits the galaxy power spectrum data down 
to small scales as well. The minimum % 2 value is 81.6 for 62 
degrees of freedom, corresponding to a probability of 4.8%. 

Fig. 1131 shows that the bias increases with increasing 
redshift, as one would expect for an old population of galax- 
ies that formed early in the first (and therefore most biased) 
overdensities. A notable exception to this trend appears to 
be redshift slice z03; however, this redshift slice corresponds 
to the region of the glitch in the luminosity-redshift distri- 
bution plotted in Fig. [5] If the median luminosity in this 
redshift slice is higher than the other slices, one would ex- 
pect a higher linear bias, consistent with what is observed. 

In order to constrain /3, we start from the definition that 
ft = f(p.j4,Q,tC)/b, where / ~ f2 m (z) 0,6 is the dimensionless 
growth factor at redshift z. Assuming that the error on / is 
larger than the variation of f2 TO with redshift, we approxi- 
mate / as a constant over the depth of the survey. We can 
then attempt to constrain / by combining all eight redshift 
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Figure 12. The measured angular power spectrum for the 8 redshift bins. The crosses show the power spectrum (and measured errors) 
of Cap I and Cap II combined, while the diamonds and stars are the measured power spectra of Cap 1 and Cap II separately. The solid 
lines are the predicted nonlinear power spectra for our fiducial cosmological model, while the dotted line shows the linear prediction. The 
dashed line is the nonlinear power spectrum for a model with negligible baryonic content. The vertical line marks the nominal nonlinear 
scale given by k = 0.1/iMpc -1 . 



slices; note that for simplicity, we ignore the correlations be- 
tween the slices and treat them as independant. The results 
are in Fig. 1141 We start by noting that the width of the 
\ 2 distribution is significantly larger than the variation in / 
with redshift, justifying our starting assumption. This is a 



direct, albeit crude, measure of Q m (z ~ 0.5) « 0.97 ± 0.53, 
consistent with our fiducial model of Sl m (z = 0.5) = 0.59. 
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Figure 14. \ 2 as a function of the dimensionless growth factor, 
/ = /3b ~ f^m: rnarginalizing over bias, for all eight redshift slices 
combined. The dashed lines mark the 1— and 2 — <r intervals, 
while the shaded region corresponds to the value of between 
z = 0.2 and z = 0.6, assuming a present day value of 0.3. The x 2 
value of 81.6, for 62 degrees of freedom has a probability of 4.8%. 
Note that we have ignored the correlations between the different 
redshift slices. 



3.6 Redshift correlations 

An important test of systematics is the cross correlation be- 
tween different redshift slices. For well separated slices, the 
cosmological correlation goes to zero on all but the largest 
scales; the detection of a correlation would imply the pres- 
ence of systematic spatial fluctuations caused by eg. stel- 
lar contamination, photometric calibration errors, incorrect 
extinction corrections etc. On the other hand, the cosmo- 
logical cross correlation is nonzero for adjacent slices due to 
overlaps in the redshift distribution, but is completely deter- 
mined theoretically by the observed auto-correlation power 
spectra and the input redshift distributions. These cross cor- 
relations therefore test the accuracy of the estimated redshift 
distributions, and in particular, the wings of these distribu- 
tions where they overlap. 

For computational convenience, we estimate the cross 
correlations with a simple pseudo-C; estimator, 



C 



21 + 1 ^ 



til ,lmO<2,lm 



(32) 



-i 



where oi,a,jm are the spherical tranforms of the galaxy den- 
sity. The pseudo-C*! power spectrum is the true power spec- 
trum convolved with the angular mask of the survey; it is 
therefore convenient to work with the cross correlation co- 
efficient, 



or- 



M * Cl 



(33) 



where M-k represents convolutions by the angular mask. 
The advantage of the cross correlation is that on scales 



Figure 17. The Fisher z-transform of the correlation coefficient 
between redshift slices separated by at least two redshift slices. 
The curves show the 1— and 3 — <r contours, given the null hy- 
pothesis of no correlations. Note that the window function is only 
approximately corrected; nearby multipoles are therefore corre- 
lated. 



smaller than the angular mask, the effect of the angu- 
lar mask approximately cancels, allowing for easy compar- 
ison with theory. It is useful to apply Fis her's z-transform 
{Kendall fc Stuardll977llPress et al]ll992h 

'1 + r N 



2 6 Vl 



(34) 



which is well described (for Z>50) by a Gaussian with mean, 



1 



loj 



(i 



1 + Ttruc 



2 VI Vtrue 

and standard deviation, 

, > 1 



+ 



Tfr 



2(N- 1) 



(35) 



(36) 



where N m (21 + 1) f s k y is the number of independent modes, 
and r tr ue is the predicted cross correlation coefficient. 

Figs. 1151 1161 1171 show the measured cross correlations 
between adjacent and more widely separated slices respec- 
tively. The absence of correlations between widely separated 
slices indicates a lack of small scale systematics common 
to the different redshift slices. The cross-correlations be- 
tween adjacent slices broadly agree with the predictions from 
the auto-correlations, although there are differences at the 
~ 10% level as seen in the plot in the lower right. There are 
two possibilities for this disagreement. The first is that vari- 
ations in the galaxy population over a redshift slice could 
cause the bias in the overlap region to differ from the value 
averaged over the entire slice. Comparing with Fig. 1131 we 
note that slice to slice bias variations of ~ 10% are consistent 
with the data. 

A second possibility is errors in the redshift distribu- 
tions. To quantify this, we model possible redshift errors by 
a shift in the median redshift. An example of this is shown 
in Fig. 1181 for the z02 and z03 slices. The figure demon- 
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Figure 15. The Fisher z-transform of the correlation coefficient between adjacent redshift slices. The diamonds and stars are the results 
for the NCC and EC respectively; the errors are 2a errors. Note that the window function is only approximately corrected; nearby bins are 
therefore correlated. The solid curve is the prediction for the cross correlation with the bias of both the slices fixed to the autocorrelation 
value, while the dotted lines show the fits allowing a variable bias. The plot in the lower right corner shows the best fit bias compared 
to the prediction from the autocorrelation. 



strates that shifting the median by 10% of the slice width 
can account for the discrepancies in the cross-correlation 
power spectrum. Furthermore, note that the corrections 
to the auto-power spectra are ~ 5%, and is principally a 
multiplicative factor that is degenerate with the bias. Fi- 



nally, the above discussion also demonstrates that the cross- 
correlation spectra are able to constrain errors in the median 
redshift at the percent level. 
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Figure 16. The same as Fig. 1151 but for redshift slices separated by one redshift bin. The overlaps at low redshifts are negligible, but 
increase at higher redshift. 



3.7 Calibration Errors 

The final systematic effect we consider is photometric cali- 
bration errors. Fluctuations in the photometric calibration 
will select slightly different populations of galaxies over the 
entire survey region, imprinting the pattern of photomet- 
ric zeropoint errors on the derived density fluctuations. One 
expects calibration errors to result in striping perpendicular 



to the drift scan direction (approximately RA). These would 
have a characteristic scale of ~ 0.22° (the width of a camera 
column), corresponding to a multipole I ~ 800, correspond- 
ing to smaller scales than those considered in this paper. 
The situation is further improved by the fact that the SDSS 
drift-scan "strips" are often broken up into several pieces 
with different photometric zeropoints, further reducing the 
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Figure 18. The change in the angular power spectrum when the 
z02 redshift distribution is shifted by Az = 0.01 [dashed], the 
z03 redshift distribution is shifted by Az = —0.01 [dot-dashed], 
and z02 and z03 are shifted by 0.005 and —0.005 respectively 
[solid] . The ratios of the angular auto power spectra are approx- 
imately 1 (±0.05), while the ratios of the cross correlations are 
approximately 1.2. 

coherence length. Thus, on the angular scales used in this 
paper, one expects calibration errors to have an approxi- 
mately white noise power spectrum. 

A useful diagnostic of photometric calibration errors 
is the cross-correlation between redshift slices with negligi- 
ble physical overlap; calibration errors will be common be- 
tween both slices. Estimating the induced cross-correlation 
requires simulations to propagate the calibration errors 
through the selection criteria and photometric redshift esti- 
mation. We simulate this by perturbing the magnitude ze- 
ropoint of each camera column and filter separately; the re- 
sulting catalogs are then input into the LRG selection and 
photometric redshift pipelines. 

Fie. ll9l shows example cross-correlations for one of these 
simulations. The lack of an observed cross-correlation argues 
for photometric calibration errors < 2%, consistent with 
other astrophysical tests of the calibration (D. P. Finkbeiner, 
private communication). The effect of such errors on the 
autocorrelation measurements is subdominant to the statis- 
tical errors. Note that the survey scanning strategy makes 
the large scale power spectrum relatively insensitive to ~ 1% 
calibration errors, the expected calibration accuracy of the 
SDSS. 



4 THE 3D POWER SPECTRUM 

Although the above power spectra are a perfectly good 
representation of the cosmological information contained in 
these data, there are advantages to compressing these eight 
2D power spectra into a single 3D power spectrum. The first 
is aesthetic; given a cosmological model, the 3D power spec- 
trum can be directly compared to theory, in contrast to the 



Figure 19. The effect of calibration errors on the cross- 
correlation power spectrum of non-adjacent redshift slices for a 
single simulation. The solid line shows the observed cross corre- 
lation (as the Fisher z-transform), while the dotted, dashed and 
dot-dashed lines show the effects of 5%, 2%, and 1% calibration 
errors. 

2D power spectra which involve convolutions by kernels de- 
termined by the redshift distributions of the galaxies (that 
contain no cosmological information by themselves). Fur- 
thermore, the 3D power spectrum directly shows the scales 
probed, and allows one to test (in a model independant man- 
ner) for features like baryonic oscillations. Finally, the 2D 
power spectrum requires computing the convolution kernels, 
making it expensive to use in cosmological parameter esti- 
mations. We however emphasize that this is (as shown be- 
low) simply a linear repackaging of the data. 

4.1 Theory 

Inverting a 2D power spectrum to recov e r the 3D 
power spectrum has been disc ussed by ISeliakJ jl998l) and 
lEisenstein fc Zaldarriagal i200lT) . An important detail where 
the two methods differ is in how they regularize the inver- 
sion. Since the 2D spectrum is the result of a convolution 
of the 3D power spectrum, it is generally not possible to re- 
construct the 3D power spectrum exactly given the 2D spec- 
trum, and one must regularize the inversion. In practice, this 
limitation is not severe, since one would normally estimate 
the power spectrum in a finite number of bands; these reg- 
ularize the inversion if the band width approximately corre- 
sponds to t he width of th e convolu tion kernel. This is the so- 
lution that lSeliakl lll99ctl presents. lEisenstein fc Zaldarriagel 
(2001) consider bands that have sub-kernel width, and reg- 
ularize the inversion by conditioning singular modes in an 
SVD decomposition. These modes are, however, given a 
large error, and so cont ain no inform ation. We adopt the 
regularization scheme of ISeliakl jl998f) . 

We start by writing the 3D power spectrum, A 2 (fc) as, 

A 2 (fc) = 6(k)Al(k) , (37) 
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where S(k) is the sum of step functions whose amplitudes are 
to be determined, while Ao(fe) is a fiducial power spectrum 
that describes the shape of the power spectrum within a bin. 
If we now describe both the 2D power spectrum, Ci, and 
and the 3D power spectrum S(k), as vectors of bandpowers, 
Eg. 1181 can be rewritten as a matrix equation, 

Ci = W<5 , (38) 

where W is the discretized convolution kernel. The solution, 
by singular valu e decomposition or normal equations (see 
iPress et alJll99l 15.4), is 

c,- 1 = W'C^W 

5 = CsW'C^Ci , (39) 

where C c ; and Cg are the covariance matrices of C; and S(k) 
respectively. 

The above discussion glossed over a number of sub- 
tleties. The first is extending this formalism for TV 2D power 
spectra. If we assume that these 2D power spectra are de- 
rived from the same 3D power spectrum, one just expands 
Ci and C c i to contain all the power spectra. However, in gen- 
eral, the 3D power spectra that corresponds to each of the 
2D power spectra could differ both in their bias and nonlin- 
ear evolution. For the latter, we divide S(k) into two sets of 
bands, linear bands with k < k n i , and nonlinear bands with 
k > k n i- We then assume that the linear bands are common 
to all TV 2D power spectra, but that there are TV copies of 
the nonlinear bands that correspond to each of the TV power 
spectra. In what follows, we assume that k n i = 0.1/iMpc -1 . 

Accounting for differences in bias over the different red- 
shift slices (as seen in Fig. 1131 is more involved. Naively 
adding TV bias parameters to Eq. 1381 destroys the linearity 
of the system. One might simply use the best fit values in 
Fig. 1131 but the fiducial model used might not correspond 
to the best fit model. We therefore use an iterative scheme 
and minimize the L 2 norm of 

Q-WS, (40) 
b 

where & is a vector of the biases (squared); these biases are 
then held constant and the inversion is performed as above. 

The next subtlety involves the choice of j3 in order to 
compute the redshift space distortions. As Fig. 1141 shows, 
these data only weakly constrain /3, and therefore we choose 
to use the linear theory prediction for /3 (more precisely for 
/ = 6/3), since the redshift space distortions only affect the 
largest (and therefore most linear) scales. 

Finally, correctly combining the different redshift slices 
requires knowing the covariance between the slices. However, 
the power spectrum estimation in Sec. I3.4l treats each slice 
independently and does not return the covariance between 
the different slices. In order to estimate the magnitude of 
this effect, we start by observing that the covariance between 
redshift slices 1 and 2 for multipole I, C(li,h) is 

Ctfi,J a )~2[£7, ia ] a (41) 

where C ( 12 is the angular cross power spectrum. Using the 
fact that the above relation is exact for full sky surveys, 
we substitute this into Eq. 1391 and use the results with and 
without these redshift correlations to scale the errors we 
obtain from the inversion. We discuss the validity of these 
approximations below. 




0.01 0.10 1.00 

k (h Mpc 1 ) 

Figure 20. The 3D power spectrum obtained by inverting the 8 
2D power spectra, normalized to the z = 0.2 power spectrum on 
linear scales (k < k n i), and uses the z = 0.2 bands for the nonlin- 
ear bands. The solid and dashed lines represent binnings Bl and 
B2 respectively, and the two power spectra are consistent. Note 
that these binnings are not independent, and must not be com- 
bined for fitting. Also shown are the nonlinear power spectrum 
using the halof it nonlinear prescription [solid, black], the linear 
power spectrum [dotted, blue], and our suggested nonlinear pre- 
scription (see below) assuming Q = 10.5 [dashed, green] for our 
fiducial cosmology. 



4.2 Results 

The result of stacking the eight 2D power spectra to obtain 
a single 3D power spectrum is shown in Fig. 1201 Note that 
the inversion process yields eight 3D power spectra that dif- 
fer on scales k > k n i = 0.1/tMpc - ; Fie. l20l shows the power 
spectrum for z = 0.2 — 0.25 slice which covers the largest 
dynamical range in wavenumber. Also note that the normal- 
ization of the power spectrum is arbitrary; we normalize it 
to the amplitude of the power spectrum at z ~ 0.2 in the fig- 
ure. Fig. 1201 shows two different binnings (hereafter Bl and 
B2) of the power spectrum interleaved with one another; the 
consistency of the estimated power spectra demonstrates an 
insensitivity to the choice of binning. 

A second assumption necessary for the inversion is the 
choice of a cosmology to convert redshifts to distances. In 
principle, the consistency between the different slices is a 
sensitive test of the cosmological model; however, the errors 
in these data are much larger than this effect. In order to test 
this, we redo the inversion with 3 different cosmological mod- 
els, and compare the results in Fig. 1211 after marginalizing 
over the bias. Note that the changes in the power spectrum 
are significantly smaller than the associated errors, while the 
errors in the power spectrum remain virtually unchanged. 
Fig. 1211 also demonstrates that the inversion process does 
not depend on the particular shape of the prior power spec- 
trum. 

Three important features of this power spectrum are: 
• Real space power spectrum : Since the individual angu- 



© 0000 RAS, MNRAS 000, 000-000 



The Clustering of Luminous Red Galaxies in the Sloan Digital Sky Survey Imaging Data 17 



n m i n 


kmax 


A 2 




6 


&8 


0.005 


0.010 


2.8639E-04 


2.2986E+00 


8.7243E-01 


0.010 


0.025 


4.4282E-03 


1.0989E+00 


1.1675E-01 


0.025 


0.040 


2.1702E-02 


8.9660E-01 


8.2658E-02 


0.040 


0.060 


5.3956E-02 


9.1448E-01 


5.8324E-02 


0.060 


0.075 


1.0630E-01 


1.0612E+00 


6.0193E-02 


0.075 


0.090 


1.5237E-01 


9.3736E-01 


6.0019E-02 


0.090 


0.130 


2.3303E-01 


1.0118E+00 


3.2957E-02 


0.130 


0.200 


4.4947E-01 


1.0281E+00 


5.4245E-02 


0.200 


0.300 


8.5115E-01 


1.2406E+00 


5.0454E-02 



Table 2. The 3D real space power spectrum (for bins Bl). The 
bands are step functions defined by k m i n < k < k m ax , the fiducial 
power spectrum by A 2 ,, and the estimated power spectrum and 
errors by 8 and as- Note that the full covariance matrix must be 
used for any detailed fitting to these data, since different data 
points are anti-correlated. 
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OS 


0.007 
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9.4449E-01 


2.8597E-01 
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6.6288E-02 


0.065 


0.080 


1.2342E-01 
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5.6484E-02 
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1.6452E-01 


9.7427E-01 


6.3003E-02 
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5.9607E-01 
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Figure 22. The ratio of the measured power spectrum to the 
linear CDM power spectrum for our fiducial cosmology (without 
baryons). As above, the solid and dashed lines represent binnings 
Bl and B2 respectively. Also shown is the same ratio for the 
nonlinear prescription, and the "no-wiggle" fit to the power spec- 
trum. The difference in \ 2 between these two models is shown for 
the two binnings. Also note the baryonic suppression of power on 
large scales, and the rise in power due to nonlinear evolution on 
small scales 



Table 3. Same as Table |3 except for bins B2. 



0.4 r , 

C °- 2 /\ /V\A ,,N . 

5: o.o ' \ / / ■-. >u — / .--^-------<- 

< . . 2 l V v 

-0.4 

2 4 6 8 10 
Bin 



1.05 : 

o 




0-95 

2 4 6 8 10 
Bin 



Figure 21. (Top) The change in the the recovered power spec- 
trum relative to the error, after marginalizing over a constant 
multiplicative bias, for different cosmologies/ prior power spec- 
trum shapes. (Bottom) The ratio of the errors relative to the 
fiducial case for the same set of cosmologies/ prior power spectra. 



lar power spectra make no use of radial information, the 3D 
power spectrum we obtain is a real space power spectrum on 
small scales, avoiding the complications of nonlinear redshift 
space distortions. Note that on length scales much larger 
than the redshift slice thickness, redshift space distortions 
cannot be neglected; however, the linear approximation dis- 
cussed in Sec. l3.1.Tl will be valid on these scales. 

• Large Scale Power: Fig. 1201 shows evidence for power 
on very large (k < 0.02ftMpc~ 1 ) scales. Marginalizing over 
bands on smaller scales, the significance of the detection 
on scales k < 0.01/iMpc -1 is ~ 2a, increasing to 5.5a for 
k < 0.02/iMpc -1 . Note that these scales start to probe 
the power spectrum at the turnover scale set by matter- 
radiation equality. 

• Baryonic Oscillations: Fig. l22l shows the 3D power spec- 
trum divided by a fiducial linear CDM power spectrum with 
zero baryonic content. The baryonic suppression of power on 
large scales, and the rise of power due to nonlinear evolution 
is clearly seen. We also see evidence for baryonic oscillations 
on small scales for both binnings, although we note that 
the power spectrum estimates are anti-correlated, making a 
visual goodness-of-fit difficult to estimate. 

To estimate the significance of these oscillations, we com- 
pare the best fit model obtained in the next section, with a 
version of the same p ower spectrum that has the baryonic os- 
cillations edited out ( Eisenstein fc Hulll998l) . The difference 
in x 2 f° r these two models suggests a detection confidence 
of ~ 2.5a or ~ 95%, assuming approximately Gaussian er- 
rors. A similar result is obtained in the next section from 
cosmological parameter fits to the baryon density. 
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Figure 23. Measurements of Qm (upper) and Q (lower) for each 
of the eight 3D power spectra. The solid lines use binning Bl, 
while the displaced dashed lines use binning B2. Note that Hm is 
insensitive to the redshift slice used, while Q depends sensitively 
on the particular choice of slice. 



5 COSMOLOGICAL PARAMETERS 

We defer a complete multi-parameter estimation of cosmo- 
logical parameters to a later paper, but discuss basic con- 
straints below. We consider a ACDM cosmological model, 
varying the matter density Qm and the baryonic frac- 
tion Qb/QM and fixing all other parameters to our fiducial 
choices. 

The principal complication to using the galaxy power 
spectrum for cosmological parameter estimation is under- 
standing the mapping from the linear matter power spec- 
trum to the nonlinear galaxy power spectrum, both due 
to the nonlinear evolution of structure and scale -dependent 
bias. We use the fitting formula proposed by ICole et ail 
(2005), 



= b 



1 + Qk 2 
1 + Ak 



(42) 



A 2 W 

where A = 1.4 is appropriate for a real-space power spec- 
trum, and b and Q are two "bias" parameters that we 
add to the cosmological parameters we estimate. Comparing 
this parametrization to a red galaxy sam ple from the Mil- 
lenium simulations l|SDringel et aljfeoOSf) . shows that this 
parametrization correctly describes the effects of scale de- 
pendent bias and nonlinear evolution up to wavenumbers 
k ~ 0.5/iMpc -1 (Volker Springel, private communication). 
We fit the data to k = 0.3/iMpc -1 . 

A second complication is that the inversion proce- 
dure of the previous section only combines wavenumbers 
< 0.1/i.Mpc -1 ; fitting the data beyond this requires choos- 
ing one of the eight redshift slices. In order to decide which 
slice to use, we estimate \ 2 on grids varying Qm, Q, and b 
for each of the eight 3D power spectra. We fix the baryonic 
density to Qt — 0.05, although allowing it to vary does not 
change the results. 
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Figure 25. The likelihood distributions for £Imi Q and Qb/QM- 
The solid line uses binning Bl, while the dashed line uses B2. As 
in Fig. 1241 we truncate at k = 0.3/iMpc — 1 . Also shown are the 
best fit values and la errors for bins Bl (top) and B2 (bottom). 



The best fit values for Qm and Q (marginalizing over 
the other parameters) , for each of the eight redshift slices are 
shown in Fig. 1231 We note that Qm and its error is insensi- 
tive to the choice of redshift slice, although Q depends on the 
particular redshift slice used. This is due to the fact that Qm 
is constrained by the location of the turnover in the power 
spectrum, and the shape of the power spectrum in the lin- 
ear regime, while Q depends on the power spectrum beyond 
0.1/iMpc -1 . In what follows, we use the redshift slice cor- 
responding to photometric redshifts between 0.45 and 0.50, 
as it corresponds to the median redshift of the full sample. 
However, we emphasize that all results below, except for the 
"nuisance" bias parameters, are insensitive to this particular 
choice. 

Fig. 1241 shows 2D projections of the (Qm, Qb/QM, Q) 
parameter likelihood space; the multiplicative bias b is 
marginalized over. The minimum x 2 values are 5.99 and 
6.94 (bins Bl and B2, respectively), for 5 and 6 degrees 
of freedom, consistent with a reduced \ 2 °f 1 P er degree 
of freedom. Bins Bl and B2 give consistent values for the 
cosmological parameters; B2 constrains all parameters (es- 
pecially Q) better than Bl because of the extra binning and 
the larger k range probed. We note that Q is correlated with 
Qm, since both QmIi and Q determine the broad shape of 
the power spectrum. An important consequence of this de- 
generacy is that an accurate estimation of Qm and its error 
requires varying Q; fixing or restricting Q will result in a 
biased Qm and an underestimation of its errors. 

Fig. 1251 shows the ID likelihoods for (Qm , Qb/QM, Q), 
marginalizing over all other parameters; the binnings are 
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Figure 24. Joint 2D likelihood distributions for fijKi Q and Cj,/f2 m , fixing h = 0.7 and n = 1, and marginalizing over the galaxy bias. 
The contours show A\ 2 = 2.3, 6.17 and 9.21. The left column panels use binning Bl, while the right column panels use B2. We truncate 
fitting at k = 0.3hMpc -1 (using the midpoint of the bins). Note that the two binnings are consistent with each other, with the B2 
binning providing slightly tighter constraints. 



again seen to be consistent. The likelihood for Qt/^M also 
allows to estimate the significance of the detection of bary- 
onic features in the power spectrum. The difference in x 2 
between the best fit model and the zero-baryon case is 5.75 
and 6.4 for bins Bl and B2 respectively, suggesting a 2.5a 
detection consistent with the model independent estimates 



made in the previous section. The signific ance of this resul t 
is similar to the results from the 2dFGRS (|Cole et al]l200ffl . 
but is weaker than the detectio n in the spectroscopic LRG 
sample feisenstein et al.ll2005bf) . 

Summarizing these results, we have 
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• For bins Bl : 

Q M = 0.27 ± 0.03 
Q b /Q. M = 0.18 ±0.03 

Q = 20.3 ±3.6 (43) 

• For bins B2 : 

Q M = 0.30 ±0.03 
Sl b /Sl M = 0.18 ±0.04 

Q = 15.6 ±2.6 (44) 

I n light of the recent WMAP results dSpergel et alJ 
2006), it is interesting to understand how the above reults 
change if we deviate from a scale-invariant primordial spec- 
trum. Minimizing \ 2 over Qui, Qt/Q]\i, and Q assuming 
n — 0.95, we find that (for bins B2), 

Qm = 0.31 ±0.03 
n b /Q M = 0.16 ±0.04 

Q = 16.3 ±2.8. (45) 

Reducing n (while keeping as fixed) boosts the power on 
large scales, but suppresses it on small scales. This results 
in a better fit on large scales, and a worse fit on small scales. 
To compensate for this, the best fit value of 0.^/Q.m de- 
creases (reducing Silk damping) while Q increases, boost- 
ing the power back up on small scales, while leaving the 
large scale power spectrum unchanged. The minimum \ 2 is 
marginally worse (7.24) than the scale invariant case. Note 
however that all the parameters are within the 1-a errors of 
those obtained assuming scale invariance. 
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Figure 26. The solid line shows the likelihood of the ratio of 
the fiducial distance to z = 0.5 to the observed distance k/ko, 
marginalizing over Q and the galaxy bias, but fixing Qi,h 2 = 
0.0223, U M h 2 = 0.127 and Q = 15.6. The dotted lines show 
slices through the 2D likelihood distribution of Q and k/ko a t 
Q = 16.5, 19.5 and 22.5 (from left to right). Note that the Q 
and k/ko are approximately orthogonal directions; varying the 
nonlinear correction doesn't change the distance scale. The best 
fit value of 0.92 ±0.06 is shown by the dashed lines. The thin solid 
line shows the k/ko likelihood for a negligible baryonic fraction; 
the distance constraint degrades to a 10% measurement. 



5.1 Distance to z = 0.5 

A potential application of the galaxy power spectrum is as 
a standard ruler. The two features of interest, the turnover 
and the baryon oscillations are determined by the physi- 
cal baryon and matter densities - Qbh 2 and Om/i 2 . Both of 
these are precisely determined by the peak structure of the 
CMB power spectrum. Therefore, in order to understand the 
sensitivity of the current measurements as standard rulers, 
we fix Qth 2 = 0.0223 and Q.mh 2 = 0.127, and vary Q and 
the comoving distance. In general, one would need to vary 
the comoving distance to each of the 8 redshift slices and 
recompute the power spectrum. However, given the S/N of 
the baryonic oscillations and turnover in these data, we sim- 
ply translate the 3D power spectrum in k with reference to 
our fiducial cosmology at the median redshift of the slice ko. 
The likelihood is in Fig. 1261 these data can constrain the dis- 
tance to z — 0.5 to 6%. Note that this is for a fixed value of 
Q,Mh 2 . Assuming a 10% uncertainty in Q,Mh 2 from current 
CMB measurements results in a ~ 2.5% uncertainty in the 
sound horizon, increasing the distance error to 6.5%. This 
must be compared to the 5% measurement of the distance 
to z ~ 0.35 measured by the spectroscopic LRG sample. 

Equally interesting is that Q and k/ko are orthogonal; 
the distance measurement does not change for different val- 
ues of the nonlinear correction. This highlights an important 
property of baryon oscillations as a distance measurement 
- it is relatively insensitive to the nonlinearity corrections 
that affect the galaxy power spectrum. 

We would also like to understand the fraction of the 



distance constraint from baryonic oscillations as opposed to 
the power spectrum shape. Fig. !26l also shows the likelihood 
for a model with a negligible baryonic fraction; the distance 
accuracy degrades to 10%, suggesting that most of the con- 
straint comes from the oscillations. 



6 DISCUSSION 

6.1 Principal Results 

We have measured the 3D clustering power spectrum of lu- 
minous red galaxies using the SDSS photometric survey. The 
principal results of this analysis are summarized below. 

• Photometric redshifts: This analysis demonstrates the 
feasibility of using multi-band imaging surveys with well cal- 
ibrated photometric redshifts as a probe of the large scale 
structure of the Universe. Accurate photometric redshifts 
are critical to being able to narrow the range of physical 
scales that correspond to the clustering on a particular an- 
gular scale, and thereby estimate the 3D power spectrum. 

• Largest cosmological volume: Using photometric red- 
shifts allowed us to construct a uniform sample of galaxies 
between redshifts z = 0.2 to 0.6. This probes a cosmological 
volume of ~ 1.5/i _3 Gpc 3 , making this the largest cosmolog- 
ical volume ever used for a galaxy clustering measurement. 
The large volume allows us to measure power on very large 
scales, yielding a ~ 2a detection of power for k < 0.01/iMpc, 
increasing in significance to ~ 5.5a for k < 0.02/iMpc. 
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• Real Space Power Spectrum: This power spectrum is 
intrinsically a real space power spectrum, and is unaffected 
by redshift space distortions on scales k > 0.01/iMpc~ . This 
obviates any need to model redshift space distortions in the 
quasi-linear regime, allowing for a more direct comparison 
to theoretical predictions. 

• Baryonic Oscillations: The 3D power spectrum shows 
evidence for baryonic oscillations at the ~ 2.5a confidence 
level, both in the shape of the 3D power spectrum, as well as 
fits of the baryonic density. We emphasize that this is only 
possible in the stacked 3D power spectrum, and therefore 
relies on accurate photometric redshift distributions. 

• Cosmological Parameters: The large volume and small 
statistical errors of these data constrain both the normal- 
ization and scale dependence of the galaxy bias. Using a 
functional form for the scale dependence of the bias moti- 
vated by N-body simulations, we fit for the matter density 
and baryonic fraction jointly, and obtain £Im — 0.30 ± 0.03 
and n h /Q, M = 0.18 ± 0.03. 

6.2 Using these results 

For cosmological parameter analyses, we recommend di- 
rectly using the 3D power spectra (binning B2), fitting 
both the galaxy bias (6) and its scale-dependence (Q) to 
k — 0.3foMpc . Electronic versions of all the power spec- 
tra, and covariance matrices used in this paper will be made 
publically available. In addition, a simple FORTRAN subrou- 
tine that returns \ 2 given an input power spectrum will also 
be made public. 



6.3 Comparison with other results 

Fig. 1271 compares the LRG power spectrum (B2 binning) , 
with the power spec trum obtained from the SDSS MAIN 
spectroscopic survey jTeemark et all2004) and the 2dFGRS 
llCole et al.ll2005f) : these three samples will be referred to as 
LRG, MAIN, and 2dF throughout this section. The solid 
and dotted lines show our nonlinear and linear fiducial power 
spectrum. Note that the normalization is arbitrary, and that 
we have not attempted to deconvolve the 2dF window func- 
tion. 

The two principal differences between these surveys 
and the data presented here is the volume probed, and 
the density of objects. As both the MAIN and 2dF are 
at low redshifts (median z ~ 0.1), the volume probed is 
< 0.05/i _3 Gpc 3 , whereas our sample probes 1.5/i _3 Gpc 3 (at 
a median redshift of z ~ 0.5) allowing us to measure the 
largest scales with smaller statistical errors, even with crude 
redshift estimates. This is clearly evident from Fig. 1271 where 
the LRG power spectrum extends to smaller k than either 
of the other two power spectra. 

On small scales, we again emphasize that the LRG 
power spectrum is naturally a real space power spectrum, 
and is unaffected by redshift space distortions. By contrast, 
the 2dF P(k) is in redshift space, and the MAIN P{k) which 
involves attempting to correct for linear redshift space dis- 
tortions. Note that the SDSS P(k) falls below the nonlinear 
power spectrum at k ~ 0.3/iMp c -1 , in line with the simu- 
lation results of iTeemark et all <l2004t) that motivated the 
discarding of k > 0.2foMpc -1 data from the cosmological 
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Figure 27. Comparison between power spectra from the SDSS 
MAIN spectroscopic survey (bottom), 2dFGRS P(k) (top), and 
the photometric SDSS LRGs (binning B2) from this work (mid- 
dle) . The normalizations are arbitrary, and have been simply cho- 
sen to separate the three power spectra. Adjacent LRG P(k) er- 
rors are anticorrelated, while the 2dFGRS errors are strongly cor- 
related. The dotted line shows the linear power spectrum for our 
fiducial cosmology, whereas the solid line is the halof it prescip- 
tion for the nonlinear power spectrum. Note that the LRG power 
spectrum fits the nonlinear power spectrum to k ~ lhMpc - . 



Survey 






SDSS MAIN 


0.297 (+0.0219, 


-0.0196) 


2dFGRS 


0.271 (+0.021, - 


-0.0187) 


SDSS LRG (Bl) 


0.260 (+0.0111, 


-0.0102) 


SDSS LRG (B2) 


0.286 (+0.0119, 


-0.0111) 



Table 4. The best fit values for Qm assuming U b h 2 = 0.024, 
h = 0.72 and a scale invariant initial perturbation spectrum. We 
use the best fit nonlinear prescription suggested by the respec- 
tive authors to fit the power spectrum to k = 0.2hMpc — 1 . The 
numbers in parentheses are the upper and lower 1 — a errors. 

parameter analysis. This is a manifestation of nonlinear red- 
shift distortions, which are particularly important given re- 
cent results that suggest that redshift distortio ns go nonlin- 
ear o n larger scales than previously anticipated dSlosar et alJ 
l200ft . 

Each of the three power spectra are consistent with the 
overall shape of the fiducial power spectrum, suggesting that 
they are consistent with each other. In order to make this 
precise and to compare statistical power, we fit for !)m and 
the galaxy bias assuming Qbh 2 = 0.024, h = 0.72 and a 
scale invariant initial perturbation spectrum. To ensure a 
fair comparison, we fit all power spectra to k — 0.2ftMpc~ 
using the best fit prescription for nonlinearity suggested by 
the authors. The results are summarized in Table[I] we find 
that all three power spectra yield consistent values for Qm- 
The LRG power spectrum, however, reduces the error by a 
factor of ~ 1.75 compared with previous results. 
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On the other hand, the SDSS LRG spectroscopic sample 
is similar to this sample. The effective spectroscopic LRG 
volume is 0.75/i~ 3 Gpc 3 at a median redshift of z ~ 0.35. 
However, the spectroscopic LRGs are sparser, with shot 
noise responsible for approximately half the statistical error 
on all scales. One can compare the S/N of the two samples 
as follows - since we are only using auto-correlations of the 
redshift slices and are ignoring correlations between differ- 
ent redshift slices, we are losing a factor ~ 7 in the number 
of modes (most of the remaining cosmological information 
is contained in adjacent redshift slices). We however gain a 
factor ~ 2 from the increased volume, and another factor 
~ 2 from the higher spatial density of objects, suggesting 
that the SDSS spectroscopic LRG sample would be a factor 
of ~ 2 greater in S/N than the photometric sample. This is 
borne out by the fact that the spectroscopic sample detects 
baryonic oscillations with A^ 2 = 11.7, while the photomet- 
ric sample has A\ 2 = 6.04, about a factor of 2 smaller. Note 
that this analysis breaks down both on the largest scales 
(where the photometric survey has more leverage because 
of the greater volume), and on scales smaller than the red- 
shift errors (where the spectroscopic sample resolves more 
modes). In principle, one could gain further by using the 
cross correlations between different redshift slices. However, 
as seen in Fig. 1181 this is very sensitive to errors in the tails 
of the redshift distribution. 

We can also compare our cosmological results with 
those obtained from the third year CMB temperature 
and polarization me asurements from the WMAP satellite 
JSpereel et all 120061. The WMAP error on Q M is domi- 
nated by the error on the Hubble constant; they obtain 
Om = 0.26^o'q3, compared with our estimate of f!ju = 
0.30 ± 0.03. They also obtain Qi,/Qm = 0.17, compared 
with Cl b /Q M = 0.18 ± 0.04. Note that WMAP favours a 
primordial scalar spectral index of n ~ 0.94; using this in- 
stead of scale invariance reduces our estimate of Vth/^M to 
0.16 ±0.04, while increasing fl M = 0.31 ±0.03. We also em- 
phasize that the errors are not directly comparable, since 
our analysis uses stricter priors. It is, however, important 
and encouraging to note that we obtain consistent results 
with a completely independent dataset. 

6.4 Future Directions 

We conclude with a discussion of the future prospects for 
photometric surveys. As of this writing, the SDSS has im- 
aged twice the area used in this paper, potentially reducing 
the errors by a factor of \f%. In addition, there are a number 
of imaging surveys planned for the near and distant future, 
the Pan-STARRS 2 and LSST 3 being two notable exam- 
ples. Both of these surveys will ultimately cover about three 
times the final SDSS area to a much greater depth, further 
increasing the volume probed. 

Baryonic oscillations are also now emerging as an im- 
portant tool to constrain the properties of dark energy. 
The tradeoff between photometric and spectroscopic ap- 
proaches to their measurement is simple - photometric sur- 
veys require wide field (> 10, OOOdeg 2 ) multi-band imaging 

2 http : //pan-starrs . if a.hawaii . edu 

3 www.lsst.org 



surveys, whereas spectroscopic surveys require large multi- 
object spectrographs. Both of these approaches are being 
actively developed, and the prudent approach would be to 
pursue both, using the results from one to inform the other. 
It is worth emphasizing that wide-field imaging surveys are 
an essential prerequisite for the other approaches (with very 
different systematic errors) to understanding dark energy, 
namely supernovae and weak lensing, suggesting a synergy 
between these techniques. 

Given the efforts underway to plan the next genera- 
tion of surveys, it is timely to compare the precision of 
the dista nce measurement w e obtain with the fitting for- 
mulae of Blake et all (1200611 . Substituting our survey pa- 
rameters into their photometric fitting formula, assuming a 
median redshift of z ~ 0.5 and a photometric redshift er- 
ror a z ~ 100/i _1 Mpc (corresponding to the redshift error at 
z ~ 0.5), we estimate a distance error of 7% as compared 
with t he actual 6% error obtained. Note that iBlake et alJ 
(2006) only use the oscillation to determine the distance, 
whereas we use the entire power spectrum. 

We can now estimate the potential sensitivity of the 
next generation of surveys. Assuming a straw-man survey 
of 20, OOOdeg 2 with a median redshift of z ~ 0.8, and pho- 
tometric redshift errors of ~ 50/i _1 Mpc, we find a factor of 
~ 5 improvement in the distance measurement, yielding a 
~ 1% measurement, the current benchmark for dark energy 
surveys. Note that this is a conservative estimate, since the 
photometric redshift accuracies assumed have already been 
achieved with the SDSS. 

In order to do so, there are a number of challenges that 
must be overcome, in addition to the brute force observa- 
tional effort required. The first is technical - this work re- 
lied heavily on having accurate, well calibrated photometric 
redshifts. Demonstrating that this is possible at higher red- 
shifts, and calibrating the redshift errors is essential. The 
second challenge is theoretical - in order to optimally use 
galaxy clustering for cosmology, we will now need to under- 
stand the connections between the physics of galaxy forma- 
tion and the observed clustering of galaxies. The hope is 
that the interplay between the two would result in a more 
complete cosmological model. 
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